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O I around a neutron star accreting at very low rates are investigated. Previous 
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The spectral properties of X-ray radiation produced in a static atmosphere 



results by Alme & Wilson (1973) are extended to the range 10 ^ < L/LEdd < 
10"'^ to include the typical luminosities, L ~ 10"^^ — 10'^^ ergss"^, expected 



H ' from isolated neutron stars accreting the interstellar medium. The emergent 

spectra show an overall hardening with respect to the blackbody at the neutron 
star effective temperature in addition to a significant excess over the Wien tail. 
The relevance of present results in connection with the observability of low- 
luminosity X-ray sources is briefly discussed. 

Subject headings: accretion, accretion disks - radiative transfer - stars: 
neutron - X-rays: stars 



1. INTRODUCTION 

Since the early seventies large theoretical efforts have been devoted to 
investigate the properties of radiation produced by accretion onto neutron 
stars, in the attempt to model the observed spectra of galactic X-ray sources. 
Up to now, the interest mainly focussed on emergent spectra for luminosities 
in the range ~ 10^^ — 10^^ ergss"^, since most of observed X-ray binaries 
have L> 10^^ ergss"^. Much fainter sources may, nevertheless, exist. It was 
suggested long ago (Ostriker, Rees, & Silk 1970), in fact, that old isolated 
neutron stars (ONSs) no longer active as pulsars, may show up as very weak, 
soft X-ray sources (L ~ 10^^ — 10"^^ ergs s~^) if they accrete the interstellar 
medium. The capabilities of present instrumentation on board X-ray satellites 
like ROSAT, make now possible the observation of such low-luminosity sources 
(e.g. Treves & Colpi 1991, Blaes & Madau 1993). 

ONSs may have low magnetic fields in which case they would accrete 
spherically. Spherical accretion onto an unmagnetized neutron star was firstly 
considered by Zel'dovich & Shakura (1969, ZS in the following) and, in more 
detail, by Alme & Wilson (1973, AW), who concentrated, however, on high- 
luminosity solutions. The lack of results for L^ 10~^ LEdd led us to reconsider 
the problem. We compute the spectrum emerging from a static spherical 
atmosphere surrounding an accreting neutron star in the case of very low rates, 
extending previous calculations by AW to L ~ 10~'^ LEdd- We found that the 
resulting spectral distributions are significantly harder than a blackbody at 
the star effective temperature. Our results may lead to reconsider the present 
estimates of the observability of ONSs, which were based on the assumption of 
a Planckian spectrum. 



2. THE MODEL 

We consider a non-rotating, unmagnetized neutron star which undergoes 
spherical accretion and is surrounded by a spherical, static atmosphere; the 
envelope material is assumed to be pure hydrogen. As the accretion flow 
penetrates into the atmosphere, protons are decelerated by Coulomb collisions 
and/or plasma interactions, their bulk kinetic energy is transferred to electrons, 
and is flnally converted into electromagnetic radiation via free-free emission. 
The input physics of our model essentially coincides with that used in previous 
studies on this subject. Equilibrium solutions were firstly provided by ZS for 
the frequency-integrated case, and by AW for the complete transfer problem. 
These investigations concentrated on the study of high-luminosity regimes 
(L> 10"^^ ergs s~^) and revealed the fundamental role played by Comptonization 
of primary bremsstrahlung photons in establishing the thermal balance of 
the emitting gas and the properties of the emerging radiation. The possible 
existence of a new "hot" solution, in which Comptonization is dominant, was 
recently discussed by TuroUa, Zampieri, Colpi & Treves (1994, Paper I). A 
detailed modelling of the interactions between the impinging flow and the static 
envelope is exceedingly complicated (see e.g. Bildsten, Shapiro & Wasserman 
1992), also because a coUisionless, standing shock can form, as originally 
suggested by Shapiro & Salpeter (1973). Following both ZS and AW (see 
also Paper I), we circumvent this problem assuming that the proton stopping 
length is a free parameter of our model, together with the total luminosity 
Loo, measured by an observer at infinity. Denoting by q the gas density, 
the column density of the atmospheric material is y = J^ gdR; the value 
yo corresponding to the proton stopping length was estimated by ZS to be in 
the range 5;^ yo^ 30gcm~^ when Coulomb collisions are dominant. As it was 
shown in Paper I, no significant expansion occurs in both "cold" and "hot" 



envelopes, so we treat the radial coordinate as a constant, equal to the neutron 
star radius R^. The heat injected per unit time and mass in the envelope is 
calculated using the expression for the energy loss rate due to Coulomb collisions 
of super-thermal protons (see Bildsten et a/.). 
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where vl = c^(l — Loo/ LEdd)/'^ is the "modified" free-fall velocity, v^ 
SkT/rUp is the proton thermal velocity, yo = (1 — 2GM^/c^R^)^/'^ is the 
gravitational redshift factor in the Schwarzschild spacetime and M* is the star 
mass. We note that, owing to the gravitational redshift, the total luminosity 
seen by a distant observer. Loo, is related to the local luminosity at the top of 
the atmosphere by L^o = 2/g-^(0). 

The run of pressure P, temperature T, monochromatic radiation energy 
density Ui, and flux F^, (both measured by the local observer) are obtained 
solving the hydrostatic and the energy balance for a completely ionized, perfect 
hydrogen gas (ZS, Paper I) coupled to the flrst two frequency-dependent 
transfer moment equations in the Eddington approximation (AW; Nobili, 
TuroUa & Zampieri 1993, hereafter NTZ) 
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Here U = f Ui, du, L = AtxRI f F^ du, LEdd = ^t^GM^c/ k^s] Kes, ^p, ^o and ki 
are the scattering, Planck, absorption and flux mean opacities, and the radiation 
temperature T^ is defined by 
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W — Wh — Wc represents that part of the injected heat, Wh, which is effectively 
converted into electromagnetic radiation within the atmosphere. Wc mimics the 
possible presence of other forms of energy transport (like convection and electron 
conduction) which are not treated in detail here (see AW for comparison); the 
actual form of Wc is discussed later on. The first two moments of the source 
function, s^ and s^, account for the exchange of energy and momentum between 
electrons and photons and, for the radiative processes we are considering, 
bremsstrahlung and electron scattering, they have the form (NTZ) 
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where B^{T) is the Planck function. The free-free opacity for a completely 
ionized hydrogen gas is 



1 _ p-hu/kT 

Kff = 1.318 X lO^^gT-^/^ g{u,T) cni^g-^ 



and a functional fit to Karzas & Latter's (1961) tables was used for the velocity- 
averaged Gaunt factor. Since equations (4) and (5) define a second order 
elliptic operator, conditions must be prescribed on the entire boundary of the 



integration domain and their form is discussed in NTZ. In particular, we assume 
that diffusion holds in the deeper layers where LTE is certainly attained and 
this automatically fixes the luminosity at the inner boundary, Li^. If all the 
energy is supplied by accretion, the total radiative flux that crosses the inner 
boundary must equal the heat transported inward by non-radiative processes, 
that is to say 
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Although the total radiation energy density and luminosity are just 
the integrals of Ui, and AiiRlFi, over frequency, we found numerically more 
convenient to derive them from the first two gray moment equations 
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Equation (9) gives trivially 
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with the condition L{yin) = Lin, while equation (10) is integrated numerically 
along with the system (2)-(5), imposing U{0) = 2F(0)/c; pressure vanishes at 
the top of the atmosphere, -P(O) = 0. 
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3. RESULTS 

Equations (2)-(5) and (10) were solved numerically by means of a finite 
diflFerences relaxation scheme (NTZ) on a logarithmic grid of 50 frequency 
bins X 100 depth zones. The adimensional frequency x = hu/kT^ was used, 
T* = T{yin), and the integration range was typically —0.7 < logx < 1.1, 
—7.6 < logy < log^i^, with yin marginally smaller than yo. A typical run 
required ~ 15 minutes of CPU time on an IBM RISK/6000. Two sets of models 
were computed, both with R^ = 12.4 km, M* = IAMq: yo = 20 gcm~^ 
and luminosities in the range 10~^ < Loo/LEdd < 0.2, yo = 5 gcm~^ and 
10~^ < Loo/ LEdd < 10~^. Our numerical method should guarantee a fractional 
accuracy better than 1% on all the variables. As a further check, the total 
luminosity, given by equation (11), was compared with the numerical integral 
of Fi, over the frequency mesh at each depth: agreement was always better than 
10%. We have also verified that our solutions with Loo^ iO~^LEdd reproduce 
almost exactly those computed by AW. For low-luminosity models, which we 
are mainly interested in, particular care must be used to handle properly the 
absorption and fiux mean opacities, since the envelope thermal balance depends 
entirely on the free-free integrated source term and the radiation spectrum 
becomes very nearly Planckian in the deeper layers. 

Results are summarized in figures 1, 2 and 3, where the emergent spectra 
and the temperature profiles are plotted for different values of Loo- In all 
these models it is yo = 20 gcm~^; solutions with yo = 5 gcm~^ show the 
same qualitative behaviour. A quite unexpected feature emerging from figures 
1 and 2 is that the spectral shape deviates more and more from a blackbody 
as Loo decreases. The model with Loo = 2.25 x 10~'^LEdd is, in fact, quite 
Planckian in shape (see also AW), showing only a moderate hard excess. On 
the contrary, solutions with Loo < 10~'^L^dd are characterized by a very 



broad maximum and by a slow decay at high energies. Comptonization is 
relatively important for L^ 1Q~'^ L Eddi similarly to what happens in X-ray 
burster atmospheres (see e.g. London, Taam, & Howard 1986). For less 
luminous models, however, non-conservative scatterings play essentially no role 
in the formation of the spectrum, as it should be expected since the temperature, 
and hence the Compton parameter, becomes lower. As can be seen from figure 
3, the temperature profile is nearly adiabatic in the inner layers where the gas is 
optically thick to true emission-absorption at all frequencies; for L^ lQ~^LEdd 
an intermediate, isothermal region is present. The sudden increase of T in the 
external layers is due to the heating produced by the incoming protons, balanced 
mainly by Compton cooling at low densities. The temperature "shock" moves 
at very low values of the column density and is nearly out of our depth mesh 
for Loo = 10"^LEdd. 

Although the hard excess present in our low-luminosity spectra may be of 
interest as far as predictions on the observability of ONSs are concerned, more 
important, in this respect, seems to be the comparison of the actual emerging 
spectrum with the blackbody at the neutron star effective temperature, 
B,y{Teff), Teff = [Loo/(47ri?^o")]^/'*, which was assumed to be the emitted 
spectrum in all previous investigations (e.g. Treves & Colpi 1991, Blaes 
& Madau 1993). It is apparent from figure 2 that model spectra with 
L^ 10~^LEdd are substantially harder than the blackbody at the star effective 
temperature. The spectral hardening can be quantified introducing a hardening 
ratio 

T 

where the radiation temperature T^ is defined in equation (6) and 
T^[Biy{Teff)] = 0.96Te//. This differs from the usual definition, 7 = Tcoi/Teff, 



where Tcoi is the color temperature, because our spectra are not always well 
fitted by a blackbody. For yo = 20 gcm~^, 7 steadily increases from ~ 1.5 
(value typical of X-ray bursters in the static phase), for L ~ 10~^ — 10~^ LEdd, 
up to ~ 2.5 for L ~ lO'^ - 10'^ Lsdd (see table 1). 

4. DISCUSSION AND CONCLUSIONS 

The significant deviation of low-luminosity spectra from a Planckian equilib- 
rium distribution could appear unexpected, since radiation is in LTE in a 
medium where the scattering depth is always much less than the absorption 
one. The source function should be Planckian and the emergent spectrum, 
formed at the thermal photosphere, should coincide with B,^{Teff). However, 
if the atmosphere develops smooth temperature and density gradients in layers 
where the medium becomes optically thin to free-free, the differential nature of 
absorption opacity plays an important role. High-frequency photons decouple 
in the deeper, hotter layers and then propagate freely to infinity, contributing to 
the high-energy part of the emergent spectral flux. At large enough frequencies, 
the observed shape of the spectrum turns out to be a superposition of planckians 
at different temperatures. This result resembles closely that of standard 
accretion disks, where the emergent spectrum shows a broad plateau due to 
the combined, thermal emission of rings at different temperatures. 

Our present result that low-luminosity spectra are harder than a blackbody 
is consistent with the previous finding by Romani (1987), who computed model 
atmospheres for cooling neutron stars. Although he considered a quite different 
physical scenario, an atmosphere in radiative energy equilibrium illuminated 
from below, the free-free opacity in his cool. He models {T^ff ~ 3 x 10^ K) 
acts much in the same way as in our faint solutions, producing a hardening 
of the spectrum. Both in Romani's and in our analysis the effects of the 
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neutron star magnetic field were ignored. An insight on the role of a strong 
B field, ~ 10^^ G, on radiative transfer was recently provided by Miller (1992), 
who showed that departures from a blackbody become less pronounced, since 
opacity is increased by magnetic effects. Finally we note that the assumption 
of a pure hydrogen chemical composition used here, is not entirely ad hoc. 
In fact, contrary to what happens in equilibrium atmospheres, such as those 
considered by Romani and Miller who allowed for different compositions, it is 
likely that metals are destroyed in the accretion flow (Bildsten et ai), leaving 
just a hydrogen envelope. 

As we already stressed, a motivation for studying the spectral properties 
of X-ray radiation coming from neutron stars accreting at low rates, stems 
from the possible detection of isolated objects fed by the interstellar gas. Their 
expected luminosities, ~ 10^^ ergss"^, could be within reach of satellites like 
Einstein and ROSAT (see Treves & Colpi; Blaes & Madau; Colpi, Campana & 
Treves 1993; Madau & Blaes 1994). Very recently Stocke et al. (1994) pointed 
out that one of the objects in the Einstein Extended Medium Sensitivity Survey 
may be actually an ONS, consistently with the original suggestion by Treves & 
Colpi. The knowledge of the emitted spectrum is fundamental in estimating the 
observability of any X-ray source and a detailed analysis of the consequences 
of our results on the detectability of ONSs with ROSAT will be presented 
in forthcoming paper (Colpi et al. 1994). Here we point out that synthetic 
spectra, being significantly harder than the blackbody at Te//, may indeed 
increase the chances of detection. In table 1 we have listed the ratios of the 
computed fiux above 0.1 keV to the blackbody one for various luminosities; the 
threshold of 0.1 keV was suggested by the sensitivity of ROSAT. The solutions 
with L = 10~^ and 10~^ LEdd can be taken as representative of the typical 
luminosities expected from ONSs embedded in the average ISM or in Giant 
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Molecular Clouds (see Colpi, Campana & Treves). As can be seen from the 
table, the ratio becomes larger than unity and the flux above 0.1 keV is from 
~ 10% to ~ 40% larger than the blackbody one for 10^^^ L^ 10"^^ ergss"^. 

Our present models could be relevant also to Soft X-ray Transients in 
quiescence, such as Aql X-1 (Verbunt et al. 1993), or in connection with low- 
luminosity globular cluster X-ray sources which emit a luminosity ~ 10~^ LEdd 
(see e.g. Hertz, Grindlay, & Bailyn 1993 and references therein). These still 
mysterious objects could be either accreting white dwarfs (e.g. cataclysmic 
variables) or neutron stars in binary systems. Their spectrum, which is still 
poorly known, could be compared with the results of our model, and deviations 
from a blackbody could be an important clue in discriminating their physical 
nature. 
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Table 1 
Characteristic Parameters for Selected Models 
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0.35 


1.01 


1.88 
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0.24 


1.02 


2.21 


2.06 


10-6 


0.15 


1.07 


2.44 


2.31 


10-^ 


0.09 


1.34 


2.64 


2.55 



^ Ratio of integrated flux to blackbody one above 0.1 keV, for models with 
yo = 20gcm-2. 

^ Hardening ratio, defined in equation (12), for models with yo = 20gcm-2. 
'^ Hardening ratio for models with yo = 5gcm-2. 
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FIGURE CAPTIONS 

Figure 1. Emergent spectra for L^o = 2.25 x 10~^ , lO"*^ , 10~^ L^^^ (full 
lines), together with the corresponding blackbody spectra at the neutron 
star effective temperature (dashed lines). 

Figure 2. Same as in figure 1 for models with L^ = 10~^ , 10~^ , 10~^ LEdd- 

Figure 3. Temperature vs. column density for different values of Loo/ LEdd- 



14 



